Theory for two-photon photoemission: transport and temperature effects 
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Using a theory which treats on equal footing transport of excited electrons and electron-phonon 
scattering, we are able to explain the temperature dependence of the relaxation time in Cu as 
recently observed by Petek, Nagano, and Ogawa. We show that the unexpected increase of the 
relaxation time with temperature results from the drastic change of the electron motion due to 
the efficiency of electron-phonon scattering: the transport is ballistic at low temperature and gets 
diffusive at room temperature. Finally, our theory also reproduces the experimental measurements 
of the two-photon photoemission (2PPE) intensity as a function of the pump-probe delay. 



Due to advances in the ultrashort laser techniques, very 
efficient tools to study the dynamics of excited electrons 
in semiconductors and metals are available. One of the 
key methods is the time-resolved two-photon photoemis- 
sion (TR 2PPE) which has the crucial advantage to al- 
low a direct measurement of the variation of the elec- 
tron distributionJKj Recent measurements on Cu have 
shown an unusual non-monotonic behavior of the relax- 
ation time as a function of energyJ3~E3 Furthermore, an 
unexpected increase-of the relaxation time with temper- 
ature was reported.til We have recently developed a the- 
ory for the dynamics of excited electrons in metals, w]aic]3. 
explicitly includes the effect of secondary electrons ,EJt-3 
in contrast— .to ab-initio calculations of the electronic 
lifetime .Ejlia It consists in solving a Boltzmann-type 
equation in the random-k approximation. Within our 
approach, we were able to reproduce the peak in the re- 
laxation time at the right energy and a linear shift with 
photon frequency, in agreement with experiments. We 
have also analyzed in detail the structure and height of 
the peak and the role of secondary electrons, especially 
focussing on the Auger contribution and the d-holc life- 
time. 

To explain the observed temperature dependence, in 
this paper we extend our model by including effects 
of transport of excited electrons out of the detection 
region and electron-phonon scattering. Note that the 
temperature dependence of the relaxation time is strik- 
ing, since on the basis of Fermi-liquid theory (FLT) one 
would expect only a very rSJj&aJl decrease of t(E, T) — 
a /[(E - E F ) 2 + 'b(k B T) 2 ]£MB Since b is of oder 1 and 
fcflT is much smaller than E — Ep, the reduction is in- 
significant for excited electrons of energy of order 1 eV. 
But if one assumes that electron-phonon collisions are 
efficient enough to change the nature of the transport 
from a regime almost ballistic at low temperature to a 
diffusive regime at room temperature, then one should 
expect a longer relaxation time at higher temperature. 
This is supported by the observation that .ballistic trans- 
port strongly reduces the relaxation timeflflO 

Let us now describe the theoretical approach we use. 
The temporal variation of the occupation of a level of 
energy E and momentum k at distance z from the surface 



is described by a Boltzmann-type equation: 
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where respectively the first term corresponds to the op- 
tical excitation, the second one describes the electron- 
electron scattering, the third is the ballistic transport 
and the fourth is the electron-phonon scattering. The 
transport term is given byEH 
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Here, v z is the z-component of the electron velocity. We 
only consider transport in the z-direction, because the 
diameter of the laser spot is much larger than the opti- 
cal penetration depth. Detailed expressions for 

in the random-k approximation (including secondary- 
electron generation) and for the optical excitation term 

are given in Ref. To derive the electron- 
opt 

phonon scattering term, we first consider the general 



expression,cilEj neglecting the transfer of energy to the 
lattice. The exchange of energy between the lattice and 
the hot electrons starts to get really effective only af- 
ter 0.5 ps,L3LJ which is not the regime we consider here. 
Also, in noble metals there is only one acoustic branch 
(1 atom/unit cell), and an upper bound for the phonon- 
energy average is Ti(uj p h) ~ ksTo, where Tjj ifj-the De- 
bye temperature. For Cu Ti{oj p h) ~ 20 meV,t3 which 
can be neglected since it is small enough compared to 
the excited-electron energy (w 1 eV). On the other hand, 
large transfer of momentum is allowed, since the Debj££ 
cut-off is of order krja w tt (a is the lattice spacing) .Ell 
This is a crucial point in our theory, because only large 
momentum transfer will efficiently change the direction of 
propagation of the electrons. Let us then express within 
these simplifications, the electron-phonon scattering con- 
tribution in the random-k approximation: 
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where T — ^-\g\ 2 (2(n) + l)p(Ep) is the electron-phonon 
scattering rate. p{Ep) is the density of states at the 
Fermi surface, f2k' denotes the solid angle, and (n) = 
{exp \h(ui p h) I '{UbT)] — I} 1 denotes the thermal average 
of the phonon occupation. Additionally we assume that 
the coupling function g is constant. Noticing that fcflT> 
h(ujph}, we get the well-known formula 



r = ^-\k B T , 



(4) 



where A = 2\g\ 2 p(E f) / (H(ujph)) i^-4h e so-called electron- 
phonon mass enhancement factor Ej The last step in eval- 
uating Eq. (|3|) is to make the substitutions f{E, k) — > 
/(-£?, v) 7 where v is the z-component of the velocity, and 
J-^[l_ /(jBjk /)]_^ Ei[ l_ /(jBjWi )]. 

The 2PPE intensity is calculated as the convolution 
of the probe laser intensity P(t) with the distribution of 
excited electrons in the vicinity of the surface x3 

/oo 
dt P(t - At) 
-OO 

dze-*/ x °°°f(E,z,t), (5) 



where A CS c = 1.6 im is the escape depth taken from over- 
layer experimentsHj The effective relaxation time r is ex- 
tracted from the 2PPE intensity as a function of pump- 
probe delay, as in experiments. 

In order to allow a direct comparison between our cal- 
culations and the available experimental data we also use 
in the calculations a laser pulse of duration 12 fs and en- 
ergy 3.1 eV. For Cu, we use 15 nm for the optical pene- 
tration depth and Vt = 1.8 nm/fs for the transport veloc- 
ity. We chose th = 35 fs for the <i-hole lifetime, since it 
was shown in a previous study that this provides a good 
order of magnitude foiithe height of the peak in the re- 
laxation time in Cu.Ej This is also in agreement with 
two independent experimental measurements which have 
suggested a lower bound of order 25 fs.EBE] We have no 
other free parameter, since the parameter A which enters 
in the electraa-jphonon scattering rate T was measured to 
be A w 0.15.BEI Note, by using Eq. M) one can already 
get an insight of the effect of the electron-phonon colli- 
sions on the transport: at T — 300 K the average time 
between two elastic collisions is r e - p h = 1/T r* 30 fs, 
in agreement with the value extracted from resistivity 
measurements {r e - p h = 27 fs).E3 However, at T = 50 K, 
we get T e - p h « 170 fs. Thus, since an excited electron 
at 1 eV has a relaxation time of order r e = 50 fs, one 
expects almost no effect on transport due to electron- 
phonon scattering at T — 50 K, because r e _ p ^ s> r e . In 
contrast, at T — 300 K, one has T e - p h < r e , and thus the 
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FIG. 1. Temperature dependence of the, relaxation time r 
(experimental data for Cu(lll) from Ref.tiJ). The dotted line 
shows the 2PPE spectrum I 2PFE (At — 0) in arbitrary units. 
The pulse duration is 12 fs and the photon energy 3.1 eV. 



phonons will be very effective in changing the nature of 
the transport. 

In Fig. [I] we have plotted both the 2PPE intensity 
I 2PPE (At = 0) and the relaxation time r at T = 300 K 
and T — 50 K as a function of energy. We compare 
our results with the experimental data of Ref. [ll]. First, 
the 2PPE intensity in both theoretical and experimen- 
tal data compare quite well. A pronounced peak due to 
transitions from the d band appears at E — Ep = 0.9 eV 
and is followed by a sharp threshold at around 1.1 eV. 
At low energy, we again observe an increase of the 2PPE 
signal. Second, the data for the relaxation time show a 
surprisingly good agreement: i) the positions of the peak 
and the dip in the relaxation time are identical; ii) the 
structure is similar, although the height of the peak is 
larger in our calculation, which could indicate that the 
d-hole lifetime Th could be smaller than the value of 35 fs 
considered here; iii) the magnitude of the change due to 
the temperature is the same. For example, at E = 0.6 eV 
and 300 K we get for both the experimental and calcu- 
lated relaxation time r exp — 175 fs and r t h — 170 fs, 
while at 50 K we find r exp = 105 fs and r t h = 110 fs. 
This agreement is surprisingly good. As expected at suf- 
ficiently high energy E — Ep > 1.5 eV the relaxation time 
is almost unaffected, since the excited-electron lifetime is 
smaller in both cases than r e _ p ^. 

In Fig. H we compare the 2PPE intensity as a function 
of the pump-probe delay with the data from Ref. [IT]. One 
can observe a strong temperature dependence, especially 
at low energy. Again we get very good agreement for 
both the quantitative and qualitative aspects. Note that 
the experimental results include a peak at At — 0, which 
is absent in the calculations because our theory does not 
include coherent effects. So one should compare the data 
for At < —20 fs. First, at low energy up to 0.7 eV, in 
both cases a clear delayed rise is observed at T = 300 K. 
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FIG. 2. Temperature dependence of the 2PPE intensity for 
different energies E — Ef as a function of pump-probe delay. 
The data on the right side are experimental results from RefJij 
The data on the left side are the theoretical results. Note the 
delayed rise at low energy for T = 300 K which leads to the 
larger relaxation time. 



Second, the magnitude of the temperature effect agrees 
surprisingly well with the experimental data. It is more 
pronounced at low energy and decreases with increasing 
energy. In a very small window of energy around 1.1 eV 
(just above the ci-band threshold), we observe a reap- 
pearance of the delayed peak accompanied by a slight 
new increase of the temperature dependence. The in- 
crease at 1.1 eV was also noted in Ref. [ll], although the 
delayed rise seems to be absent. It is also interesting to 
remark that such a rise could be observed for the first 
time, since a sufficiently short laser pulse was used. This 
rise is in fact a signature of the presence of secondary 
electrons. 

To illustrate the reason why the effect of temperature 
is so strong, let us analyze the motion of the excited- 
electron distribution as a function of time. As a measure 
for the penetration of electrons into the bulk, we define 
the average distance from the surface, (z(t)): 
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where N(z, t) = / °° dEp(E)f(E, z, t) is the average num- 
ber of excited electrons at distance z and time t. In Fig. |^, 
we have plotted (z(t)) at different temperatures, T = 0, 
50, and 300 K. Clearly we observe that at T = and 
50 K, the motion of the excited electrons is ballistic. At 
T = 50 K, a small deviation from the linear behavior 
appears around t ~ 200 fs, which is of the order of mag- 
nitude of T e - P h = 170 fs. We get for the average velocity 
A(z)/At w Vt/2, where the factor 1/2 can be understood 
easily by considering the average of the velocity in z direc- 
tion. However, the nature of the motion has drastically 
changed at T = 300 K: the motion is now diffusive. We 
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FIG. 3. Average distance (z) of excited electrons from the 
surface after laser excitation. Note that at T = and 50 K 
the transport is ballistic, while at T = 300 K it is diffusive 
due to elastic electron-phonon collisions. The dashed line is 
a fit using (z) = WDtt — to). 



illustrate this by fitting the data with (z) = y/D(t — to), 
which is expected in the case of diffusive motion. An 
offset to is introduced in order to take into account the 
finite duration of the laser pulse generating excited elec- 
trons. We get D — 32 nm 2 /fs, which agrees very well 
with the expression for the electronic diffusion coefficient 
D = Vf/ e /3 = 29 nm 2 /fs, where l e — VtT e - p h is the elec- 
tronic mean free path. It is interesting to note that at 
t = 0.5 ps the excited electrons have already reached 
an average distance of 120 nm, about ten times larger 
than the optical penetration depth. This is in agreement 
with a value of 100 nm used to describe the initial spatial 
distribution of excited electrons in the two-temperature 
model.Ej Such a model does not describe the thermaliza- 
tion of the electron gas and starts to be valid only after 
t > 0.5 ps. Thus, one has to add (z(t = 0.5 ps)) to the 
optical penetration depth. 

To conclude, we have presented a theoretical model 
including both transport and electron-phonon scattering 
which is able to reproduce the temperature dependence 
of the relaxation time. It is shown that this variation 
is due to a drastic change of the excited-electron motion, 
which is ballistic at low temperature and gets diffusive at 
room temperature. Note that the correct order of magni- 
tude of the temperature effect is obtained without using 
free parameters: the electron-phonon scattering rate was 
directly taken from experimental data. Furthermore, we 
also observe at low energy a delayed rise in the 2PPE 
intensity and an increase due to temperature in agree- 
ment with the experimental results. As a final remark, 
we have provided a very efficient method to describe 
the excited-electron dynamics in the short-time regime, 
where electron-phonon energy transfer is negligible (be- 
low w 0.5 ps). The results obtained here indicate that 
it is promising to use our extended model to study other 
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problems, in particular the dynamics of excited electrons 
involving transport effects in the regime where electron- 
lattice energy transfer becomes relevant. So far, the only 
method available for this purpose is the two-temperature 
model, which is not always reliable.E-3 The present the- 
ory is also suitable for the study of thin films, which offer 
the possibility to probe the effect of transport. Our pa- 
per shows how new information about the dynamics of 
excited electrons in solids can be extracted from 2PPE. 
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